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Abstract 

We analyze scheduling algorithms for multiuser communication systems with users having multiple 
antennas and linear receivers. When there is no feedback of channel information, we consider a common 
round robin scheduUng algorithm, and derive new exact and high signal-to-noise ratio (SNR) maximum 
sum-rate results for the maximum ratio combining (MRC) and minimum mean squared error (MMSE) 
receivers. We also present new analysis of MRC, zero forcing (ZF) and MMSE receivers in the low SNR 
regime. When there are limited feedback capabiUties in the system, we consider a common practical 
scheduling scheme based on signal-to-interference-and-noise ratio (SINR) feedback at the transmitter. 
We derive new accurate approximations for the maximum sum-rate, for the cases of MRC, ZF and 
MMSE receivers. We also derive maximum sum-rate scaling laws, which reveal that the maximum 
sum-rate of all three linear receivers converge to the same value for a large number of users, but at 
different rates. 
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1. Introduction 

The use of multiple-input multiple-output (MIMO) antenna systems in multiuser environments 
has recently gained considerable attention. Much progress has been made on characterizing 
the fundamental limits and developing signal processing schemes for both the multiple-access 
channel and broadcast channel scenarios (see, eg. [1,2]). In this paper, we consider the broadcast 
channel, which embraces a large number of important application scenarios; for example, a base- 
station serving multiple users in a cellular environment, or a relay node serving multiple network 
nodes in an ad-hoc network. 

For systems with perfect feedback capabilities, i.e. when complete short-term channel state 
information (CSI) of all of the users are available to the transmitter, it is well-known that the 
optimum solution for achieving the capacity is to employ dirty paper coding (DPC). How- 
ever, while DPC is capacity achieving, it also has prohibitive complexity. This, coupled with 
the stringent and potentially overwhelming feedback requirements, render DPC unsuitable for 
practical implementation. Various low complexity alternatives have been employed and studied 
extensively (see, eg. [1,3,4]), however these schemes all require perfect knowledge of the CSI 
at the transmitter. For increasing number of users and antennas, the amount of CSI feedback 
becomes large, and these schemes become impractical. In addition, inaccurate CSI feedback may 
lead to significant performance losses [5]. 

In this paper, we are interested in the practical scenarios when there is limited or no feedback of 
CSI at the transmitter. In this case, efficient scheduling algorithms are required for coordinating 
transmissions to the different users. Appropriate scheduling approaches depend on the feedback 
capabilities between the receivers and the transmitter, and their performance can be measured 
in terms of the maximum sum-rate. 

For systems with no feedback capabilities, an obvious and practical scheduling approach is to 
employ a simple round-robin scheduling algorithm; i.e. a time-division multiple access (TDMA) 
scheme. The maximum sum-rate in this case becomes equivalent to the maximum sum-rate of 
a single user MIMO system. The single user MIMO maximum sum-rate has been extensively 
studied for optimal receivers (see, eg. [6,7]), and zero forcing (ZF) linear receivers (see, eg. [8, 
9]). Other results include the particular case where the number of transmit antennas is less than 
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or equal to the number of receive antennas in the high signal-to-noise ratio (SNR) regime, for 
which the maximum sum-rate has been derived for ZF and minimum mean squared error (MMSE) 
receivers in [10, Eq. 8.54]. In this paper, we derive new exact maximum sum-rate expressions for 
both the maximum ratio combining (MRC) and MMSE receivers. We also derive new maximum 
sum-rate expressions in the high SNR regime for MMSE receivers in the case where the number 
of transmit antennas is greater than the number of receive antennas, and for MRC receivers in 
the case of arbitrary antenna configurations. We then investigate the maximum sum-rate in the 
low SNR regime for each of the three linear receivers. 

When limited feedback is available, the scheduling algorithm can be more sophisticated. In 
this paper we consider the common case of systems which feed back signal-to-interference-and- 
noise ratios (SINRs) and use opportunistic scheduling. In particular, we examine systems where 
the SINR (for each user) can be measured from each transmit antenna separately, and then fed 
back. In the scheduling algorithm each transmit antenna sends an independent data stream to the 
user with the largest corresponding SINR. Of particular interest is systems employing practical 
linear receivers. For the specific case of ZF, previous related results have been presented in 
[11], where an upper bound to the asymptotic maximum sum-rate scaling law was derived, and 
in [8], where an exact (albeit complicated) non-asymptotic expression was derived. Moreover, 
for the specific case of four transmit and two receive antennas, [12] presented maximum sum- 
rate scaling laws for the MRC and MMSE receivers. In this paper, we generalize these prior 
results, by considering the maximum sum-rate for arbitrary antenna configurations, and deriving 
new maximum sum-rate results for MRC, ZF and MMSE receivers. We present new accurate 
approximations for the maximum sum-rate, and derive exact maximum sum-rate scaling laws. 
We show that all three linear receivers converge to the same asymptotic maximum sum-rate as 
the number of users grows large, however the speed of convergence is different. Our results are 
confirmed through comparison with Monte Carlo simulations. 

Throughout this paper, we denote E[-] as expectation, {-Y as conjugate transpose, Tr(-) as 
matrix trace, Ijv as a x identity matrix and O^xm as a x M zero matrix. Also, ® 
denotes Kronecker product and CA/'j\/,Ar(M, V) represents an M x A^ matrix- variate complex 
Gaussian distribution with M x N mean matrix M and MN x MN covariance matrix V. 
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II. System Model 

We consider a multiuser MIMO broadcast channel scenario, where a transmitter has Nt 
antennas while each of the K users has Nr antennas. Each antenna at the transmitter is used to 
independently transmit a stream of data to a particular user, as determined by the two scheduling 
algorithms described in the next two sections. The received vector for the zth user can be written 

as 



HjX + rij 



(1) 



where Hj ~ CN'Nr,Nt {^NrxNti^Nr ® Ia^J is the Nr x A^^ Rayleigh fading channel matrix from 
the transmitter to the 2th user, x is the Nt x 1 transmit symbol vector with Efx^^x] = P, and 
rij ~ CMn,.,! (Oat^xi, Nol^r) is the additive white Gaussian noise vector. If we assume that the 
data stream for the ith user comes from the A;th transmit antenna, we can write the received 
vector in ([T]) as 

Nt 

Yi^k = iii,kXk + \jXj + Ui (2) 

where hij is the jth column vector of Hj and Xj is the symbol sent from the jth transmit antenna. 
We assume that user i perfectly estimates their own channel matrix, Hj. We also assume that 
the total power budget is distributed equally across the different transmit antennas, such that 
Ed^j p] = for j = 1, . . . , Nt. To recover the desired symbol, we multiply the received signal 
by a 1 X A'r received weight vector w], which for MRC, ZF, and MMSE linear receivers, and 
their resulting received SINR, are given respectively by 



i,k 



7i,A:,MRC 



ZF 



-2- V^* Ih"^ h- -12 + Ih'^ h- 1.1 

P 



Nt 



and 



'^IlMSE 



7i,fc,MMSE 
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(3) 
(4) 



Nt 



KKt^ 

Nt 



l-Nr 



H,k 



(5) 
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where p = is the average SNR, g^. is the kth row of (H|Hi)"^Hj, [Z]k^k is the (A;, A;)th 
element of Z and K = Hj^*^^ is the A^,. x (A^^ — 1) matrix with the same elements as Hj but with 
the kth column removed. Note that throughout this paper, we require A^,. > A^^ when dealing 
with ZF receivers, whilst we consider arbitrary antenna configurations for MRC and MMSE. 
It is convenient to define the normalized SINR as follows 

A 



X 



P 



(6) 



where 7 is the instantaneous SINR given for MRC, ZF and MMSE in ([3]), dH) and © respectively. 
For the ZF and MMSE receivers, the cumulative distribution function (c.d.f.) of X for an arbitrary 
user is given by [8] 



Nr-Nt k 
-X \^ ^ 

^ k\ 

k=0 



and [13] 



Nr-l 



Nt 



X 



Nt-1 ^ f^''^ 



(V) 



(8) 



k=0 



respectively, where 



Pk 



_P_ 

Nt 



\ k-p J 



(9) 



p=max(0,fc-Art+l) P- yjTt 

For MRC, the c.d.f. of X has been previously derived in [14]. However the resulting expression 
is complicated, involving multiple summations over sets, and is therefore not easily amenable to 
further analysis. In Appendix [Al we show that the SINR c.d.f. admits a simpler representation 
given by 



1 



Nt 



Nr-l k 
k=0 p=0 



1 + ^ X 



(10) 



where 



(^p,k 



(Nt+p~2\ ( _p_ 
\ p > \Nt 



(11) 



{k-p)\ 

The corresponding probability density functions (p.d.f.) for ZF, MMSE, and MRC are obtained 
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by taking the derivative of dV]), dH]), and (flOl) respectively, and are given by 



Q-X^Nr-Nt 



and 



= , (12) 

f 1 + j^xj k=o V V / / 

Nr~l k f._i . . X , 



(14) 

In the next two sections, we analyze the maximum sum-rate of the three linear receivers using 
two different scheduling algorithms which differ depending on whether or not there is SINR 
feedback from the different users to the transmitter. We present exact and approximate results, 
and compare the performance of the three linear receivers. Note that we assume all users employ 
the same receiver structure (ie. ZF, MMSE, or MRC). 

III. Scheduling Without Feedback 

For multiuser MIMO systems with no feedback capabilities, we consider a simple scheduling 
algorithm where each transmit antenna is assigned to users in a round robin manner. This includes 
any scenario where each user is served by each transmit antenna for equal time periods and the 
data streams received by the users are decoded independently. We consider the practical scenario 
where the number of users is greater than the number of transmit antennas. The maximum sum- 
rate under this scheduling scheme is given by 



log. [l + ^X 



(15) 



where the distribution of X is given in (|7]), ^ and (flOl) for the ZF, MMSE and MRC receivers 
respectively. Note that this is statistically equivalent to the maximum sum-rate of a single user 
MIMO spatial multiplexing system with linear receivers. Although these systems have been 
studied extensively, there still remains important open work in this area. Specifically, in this 
section we present new exact and high SNR closed-form expressions for the maximum sum- 
rate in (fT5l) with MRC and MMSE receivers. We also present a new analysis of the maximum 



sum-rate for each of the three linear receivers in the low SNR regime. 

A. Exact Maximum Sum-Rate Results 

The exact maximum sum-rate of MIMO systems with ZF receivers has been derived previously 
(see eg. [8,9]), and is given by 

/ AT \ Nr-Nt + 1 Nr-Nt + 1 / \ " / ]\J \ 

R^iM=^og,^[-j) N,el^ E [irj ^[n-Nr + N,-l,-^j (16) 

where r(-, ■) is the incomplete gamma function [15]. For MMSE and MRC receivers however, 
the corresponding expressions are not available. To derive these new results, we find it useful 
to re-express the maximum sum-rate expression in ([T?]) . as given by the following lemma. 
Lemma 1: The round robin maximum sum-rate can be written as 

R"ip) = p\og, 2 r ^-^^^dx . (17) 
Jo 1 + 

Proof: The proof follows by applying integration by parts to ([T5] ). ■ 
In Appendix |Bl we employ Lemma [D to derive the exact maximum sum-rate of MRC receivers 
as follows 

/ AT \ Nt-1 Nr-1 k 

RmkM = loge 2 [jj N.e"^ ^ 

^s(:)<-''-'(#M-— 'f)- - 

Using a similar approach (we omit the specific details), it can be shown that the exact maximum 
sum-rate of MMSE receivers is given by 

/ AT M'-l 

^mmse(p) = loge 2 [-jj N.e"^ ^ 

Fig. [T] compares the analytical maximum sum-rate for MRC, MMSE and ZF based on (fTSl) , 
([T9l ) and [8, Eq. 6] respectively. Monte Carlo simulated curves are also presented for further 
verification. As expected, we see that the maximum sum-rate of MMSE is greater than for MRC 
and ZF for all SNR values. We also see an SNR crossover point, such that the maximum sum- 
rate of MRC is greater than ZF below the crossover point and less than ZF above the crossover 
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point. This crossover point occurs because the SNR of MRC and ZF approaches that of MMSE, 
the optimal linear receiver, at low and high SNR respectively. 

We numerically calculate this ZF-MRC SNR crossover point in Fig. |2l for different antenna 
configurations using the maximum sum-rate expressions for MRC in (fTSl ) and ZF in [8, Eq. 6]. 
We see that for a fixed Nr, increasing A^^ increases the crossover point, however for a fixed 
Nt, increasing A^^ decreases the crossover point. We also see that the performance of MRC is 
better than ZF for a number of practical scenarios. For example, when Nr = Nt = 4, the MRC 
maximum sum-rate is greater than the ZF maximum sum-rate for SNR values as high as 8 dB. 

We now compare the maximum sum-rate of MRC and MMSE using the analytical expressions 
in (fTSi) and ( fT9l ) respectively. Figs.[3l|4]and[5]show the percentage of maximum sum-rate achieved 

i?''"' (p) 

by using an MRC receiver with respect to the MMSE receiver (i.e. mrcvp^ ^ j^qq-j different 



antenna configurations. We see in Fig. |3] that when A'^,. < Nt, the maximum sum-rate of MRC 
is comparable to MMSE, and approaches the MMSE maximum sum-rate for increasing Nt. For 
example, for Nj. = 4 and Nt = 7, the MRC maximum sum-rate achieves over 80% of the MMSE 
maximum sum-rate for SNR ranges as high as 5 dB. On the other hand, when Nr > Nt, we 
notice from Fig. |4] and Fig. |5] that increasing the number of transmit antennas has the detrimental 
effect of reducing the percentage of maximum sum-rate achievable. This can be explained by the 
fact that for the MMSE receiver, when A^^. > A^t, there are enough receive antennas to cancel the 
interference from the transmit antennas, hence adding more transmit antennas does not decrease 
the SNR as much as it does for the MRC receiver. 

B. High SNR Analysis 

At high SNR, for antenna configurations with A^^ > Nt, expressions for the maximum sum-rate 
of MMSE receivers have been derived in [10, Eq. 8.54]. For Nr < Nt however, corresponding 
results are not available. In this case, the maximum sum-rate of MMSE receivers approach a 
constant at high SNR. In Appendix O we derive this constant as follows' 



-MMSE 




(20) 



f(x)=o(g(x)) means that lim^j^oo = 



0. 
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Considering MRC receivers, for A^^ > 1 the maximum sum-rate also approaches a constant at 
high SNR. Following a similar approach to that used for the MMSE receiver, we evaluate this 
constant as follows 

For Nt < Nr, the maximum sum-rate of ZF and MMSE scales logarithmically with the SNR 
[10, Eq. 8.54] instead of approaching a constant, and hence always performs better than MRC as 
expected. From (|20|) and (|2T1) . this can also be shown to be the case when Nt > Nr, where the 
maximum sum-rate of MMSE outperforms MRC. However, (l20l) and (|2TI) reveal the interesting 
fact that for high SNR, if N^ is kept fixed and Nt oo, the maximum sum-rate of MRC 
approaches that of MMSE. This can be intuitively explained by noting that for MMSE, when 
Nt > Nr, there is not enough degrees of freedom for the receiver to cancel the interference 
imposed by the transmit antennas. Therefore, as A^^ grows large, both MMSE and MRC suffer 
significant interference penalty. As such, for sufficiently large Nt, the interference cancelation 
capabilities of MMSE become insignificant, and the performance of the two receivers coincide. 

The accuracy of (l20l) and (|2T1) is confirmed in the "Analytical (High SNR)" curves of Fig. 
|6l where we plot these expressions (without the o(-) terms) as well as Monte Carlo simulated 
curves. For all cases considered, we see that the Monte Carlo simulated curves converge to 
the analytical asymptotic results in the high SNR regime. In addition, we also plot the exact 
"Analytical" curves from (fTSi) and (fT9l) . and see an exact match with the Monte Carlo simulated 
curves. 

C. Low SNR Analysis 

We now investigate the maximum sum-rate of the three linear receivers in the low SNR regime. 
The simplest approach is to derive a first order Taylor expansion of R{p) as p ^ 0. It is shown 
in [16] however, that such an approach does not adequately reveal the impact of the channel, 
and may in fact lead to misguided conclusions. Moreover, as discussed in detail in [16], in the 
low SNR (or wideband) regime it is more appropriate to investigate the maximum sum-rate in 
terms of the normalized transmit energy per information bit, jr, rather than SNR. In this case 
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the maximum sum-rate representation is given for low ^ levels by the following expression 



^'^)-^°" 3dB 











No 


dB 


^Omin 


dB 







Eb 




dB 


^ ' U mm 



(22) 



dB/ 

where ^ . is the minimum normalized energy per information bit required to convey any 
positive rate reliably, and 5*0 is the wideband slope in bits/sec/Hz/(3 dB) at the point ^ . . 
Note that and -1^ can also be used to determine the bandwidth of a system designed to 
achieve a certain rate subject to power limitations [16]. 

Now as p ^ 0, it can be shown from ^ and ^ that the SINR, and hence the maximum 
sum-rate, of the MMSE and MRC receivers coincide. In Appendix |D] we derive closed-form 
solutions for and 5*0 for both MRC and MMSE as 

No min 

^^MRC/MMSE jj^2 
^Omin 

and 



(23) 



crMRC/MMSE _ 

° 2Nt + Nr 

respectively. For ZF receivers, we derive . and Sq as 



(24) 



mm 



In 2 Eb^^ In 2 

< — = < In 2 (25) 

Nr - A^Omin iV, - iV^ + 1 " 



1 < ^z.^2iW^^^ ^ 2NJ^ 

- ° Nr-Nt + 2 - Nr + 1 



and 



respectively. The proof is omitted due to space limitations, but involves solving (I5TI) using the 
p.d.f. in (fT2l) . along with some algebraic manipulation. 



We see from m and m that 1^^^^/^^^'^ < , and from m and m that ^^RC/mmse > 

Jvomin — ^Omin — 

Sq^; in both cases, with equality if A^^ = 1. Therefore, as expected, whenever A^^^ > 1, the 
low SNR maximum sum-rate of the MRC and MMSE receivers always exceeds the maximum 
sum-rate of ZF. The new results given by (|23l)-(|26l) precisely quantify this maximum sum-rate 
difference. In addition, it is interesting to note that ^^^'-'/^^^^ depends only on the number of 

ZP 

receive antennas, also observed in [17] for optimal receivers, whereas jl^^.^^ depends on both the 
number of transmit and receive antennas. In fact, we see that increasing the number of transmit 



antennas for a fixed number of receive antennas has the negative impact of increasing 
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^Omin' 



IV. Scheduling With SINR Feedback 



For multiuser MIMO systems with feedback capabilities, we consider a low-rate feedback 
approach where each user sends to the transmitter a set of Nt instantaneous SINR values. 
These correspond to the instantaneous SINRs at the output of the linear receiver; one for each 
transmit antenna. Specifically, depending on the particular linear receiver employed, the SINRs 
are calculated based on the average SNR and the channel matrix, according to either ([3]), dH) or 
([5]), for MRC, ZF, and MMSE respectively; assuming the signal from one antenna is the desired 
signal while the signals from the other antennas are interference. 

We employ an opportunistic scheduling algorithm which selects, for each transmit antenna, 
the user with the maximum SINR. Note that each user can be served by more than one transmit 
antenna, ie. multiple data streams can be sent to a single user. The maximum sum-rate can be 
written as 



where X^ax = maxi<j<;^ Xi and fx and Fx are given for MRC, ZF and MMSE in (l7])-(fT4l). 
In the following, we derive new accurate approximations for (ITTI) for MRC, ZF and MMSE 
receivers. Note that, for the case of ZF, exact expressions have been derived previously in [8]. 
Those results, however, incurred high computational complexity which increased significantly 
with the number of users and, as such, would take time comparable to Monte Carlo simulations 
to produce meaningful plots. As such, we are motivated to consider deriving simpler, more 
computationally-efficient, accurate approximations for the ZF receiver as well. 

A. Maximum Sum-Rate Approximations 

To obtain closed form expressions for the maximum sum-rate of the linear receivers, we can 
directly evaluate the integral in (l27l) by first substituting the c.d.f. expressions in (|7]), dS]) and 
(flOl) . along with the corresponding p.d.f. expressions in (fT2)) - (fT4l) into (l27l) . We can then apply 
the multinomial theorem and evaluate the resulting integral using standard identities from [18]. 




(27) 
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However, although this approach produces an exact closed-form solution, the resulting expression 
has a complexity which increases significantly with the number of users K, and would take time 
comparable to Monte Carlo simulations to produce meaningful plots for even small K values (eg. 
K = 10). As such, we are well motivated to consider accurate approximations for the maximum 
sum-rate, which have much lower computational complexity requirements. Our approach is to 
approximate (|27l) by the following 

i?--(p) ^ N, log, + ^ g{K)\ (28) 



where^ 

giK) = F^' (- (29) 

with the distribution of X given in ©, dH]) and ([TO]) for the ZF, MMSE and MRC receivers 
respectively. This expression is obtained by first noting that for a random variable Y following a 
symmetric distribution, with c.d.f. Fy(-) and largest order statistic Yk, the following inequality 
holds [19] 

HYk] < Fy' ( . (30) 

To obtain (l28l) . we substitute Y = Nt log, ^1 + ^vT^) ^^^'^ (l30l) . and apply simple algebraic 
manipulation. 

We note that this general upper bounding approach has been shown to be tight for a Gaussian 
distribution [20]. We would therefore expect it to be tight in (|28l) also, at least at low SNR 
values, upon noting that for low SNR we have log2 ^1 + ^fr^) ~ Nt'^' which has been shown 
to converge to a Gaussian distribution as Nt, Nr ^ oo for MRC, ZF and MMSE [21,22]. This 
expected tightness is confirmed in our numerical results. 

We will focus on evaluating (l28l) . and show that it yields an accurate approximation for (ITtI) . To 
calculate (l28l) . we require the inverse function of Fx{-) in g{K). Unfortunately, for ZF, MMSE 
and MRC, the c.d.f.s in ([7]), dS]) and (flOl) respectively, are in a form for which the inverse c.d.f. 
is difficult to obtain. We thus employ a simple numerical algorithm, similar that considered in 

^Note that throughout this paper, we will use F^^ to denote the inverse function of F, and to indicate the A^th power. 
There is no contradiction in this notation since we only consider positive powers in this paper. 
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[23], to evaluate g{K) = F^^ I — _J:k-i i I • In particular, we start with go{K) = 1, and iterate 
the following expression 

g,^,{K) = In {N{K)) + In (1 - Fx {gj{K))) + g,{K) (31) 

until convergence, where N{K) = 1 + e^^=i^^, and Fx{-) is given in dV]), ([8]) and (flOl) for 
ZF, MMSE and MRC receivers respectively. Note that the algorithm will eventually converge, 
because additional iterations of (|3TI) results in additional nested logarithms, which contribute 
less and less to the outermost logarithm, ie. the logarithm at the first iteration at j = 0. 

Simulations indicate that (|3TI) has a fast rate of convergence. The exact rate, however, depends 
on the number of antennas and linear receiver used in the system. Table H shows the absolute 
error percentage versus the number of iterations, for the MMSE receiver and for different antenna 
configurations. The absolute error percentages were calculated by using the algorithm in (|3TI) 
and numerical results from the exact c.d.f. in ([8]). We see that the convergence rate is fast in all 
cases. 

We denote gcom{K) to be the converged value from (|3TI) . Substituting this into (l28l) in place 
of g{K) gives an approximation to maximum sum-rate as follows 

R^^^p) ^ Nt log2 (l + |- • (32) 

Fig. |7] compares the analytical approximations for MRC, ZF and MMSE with Monte Carlo 
simulated curves. The analytical approximation curves were obtained using (|7]), ([8]), (fTOl) , (|3T|) 
and (|32|) . We see that these curves accurately approximate the Monte Carlo simulated curves in 
all cases. 

B. Large User Analysis 

We now examine the maximum sum-rate for MRC, ZF and MMSE, in the large-K regime . 
Theorem 1: For sufficiently large K, the maximum sum-rate for the MRC, ZF, and MMSE 
receivers is given by 



= N, log2 In K + N, logj p - N, log^ Af, 



AT. log, ^H^^.^^L^PiN.-N.mnlnKY 
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where a is given for MRC, ZF and MMSE as follows 



( (Nt\Nt-l 



MRC 



= < 



1 



ZF 



(34) 



a 



{Nr-my. 

Nr-l 



JVr-l-p. 



E 



MMSE . 



Proof: See Appendix El 
Corollary 1: For even larger K, the maximum sum-rate (1331 ) becomes 



^p=ma.x{0,Nr~Nt) 




(35) 



Note that for the case of ZF receivers, an upper bound to (l35l) was derived in [11], and for 
the case of MRC and MMSE receivers with A^i = 4 and Nr = 2, the exact maximum sum-rate 
scaling law was reported in [12]. As such, our expressions generalize these prior results by 
giving exact scaling laws for all three receivers, which apply for arbitrary numbers of antennas. 
Our results also yield important insights into the fundamental differences in terms of maximum 
sum-rate scaling for each of the three linear receivers. 

From Corollary [H we see that the maximum sum-rate for all three receivers follow the same 
asymptotic scaling law with respect to the number of users. The key difference, however, is 
demonstrated by Theorem [H where the rate of maximum sum-rate increase with K is seen to 
be a decreasing function of a. That is, although the linear receivers follow the same asymptotic 
scaling law, the speed at which that law is approached varies with a, which depends on the 
particular receiver structure. 

In general terms, we see from (1331) that for large but finite numbers of users K, the maximum 
sum-rate is a monotonically increasing function of a. From (|34l) . we clearly see that cimmse > 
ttMRc and ommse > ctzF- We also have auKC > "zf if 



and azF >= omrc otherwise. Therefore, for large K, (|36l) represents the SNR crossover point 
between the maximum sum-rate of ZF and MRC receivers, discussed in Section ITlI-AI We note 
that the right-hand-side of (|36l) is a decreasing function of the number of receive antennas. 
This is consistent with the behavior of the corresponding ZF-MRC SNR crossing point for the 




(36) 
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round-robin scheduler, observed in Section IIII-A[ 

Note that although this convergence behavior of the three linear receivers is not readily apparent 
from the results in Fig. Ul for the finite numbers of users considered, it can be clearly seen in Fig. 
[8l where we plot the difference in maximum sum-rate between the MMSE and MRC receivers 
(i.e. -Rmmse(p) ~ -^MRc(p)) ^ function of the number of users using (l32l) . We clearly see that 
the difference in maximum sum-rate is a decreasing function of the number of users. 

The fact that all three linear receivers converge in the limit as A' ^ oo can be intuitively 
explained as follows. The performance in each case is limited by the interference from the 
other transmit antennas. As K oo, the probability of the opportunistic scheduling algorithm 
choosing a user such that the interference is negligible approaches one. In this interference-free 
scenario, all linear receivers thus have the same performance, and the SINRs all converge to the 
same value, i.e. 7mrc ^ 7zf ^ 7mmse as K oo. 

We see in Fig. |7]that the opportunistic scheduler can provide a significant increase in maximum 
sum-rate compared to a round robin scheduler (i.e. corresponding to the point K = 1); at the 
expense of requiring more SINK feedback. It is also worth noting that when the users are 
stationary, fairness may become an issue for the SINR-based scheduler, since the transmitter is 
likely to send to a subset of users for the majority of the time. When the users are sufficiently 
mobile, however, such that the channel realizations change for each transmission period, then on 
average each user will be served for approximately equal amounts of time under the opportunistic 
SINR-based scheduling approach. 

V. Conclusion 

We have analyzed the maximum sum-rate of MRC, ZF and MMSE receivers in MIMO 
multiuser systems. We considered two scheduling algorithms, which differ depending on whether 
or not there is SINR feedback at the transmitter. When there is no feedback, we considered a 
simple round-robin scheduling algorithm, deriving new exact closed-form maximum sum-rate 
expressions for the MRC and MMSE receivers for arbitrary SNRs, and simplified results in the 
high SNR regime. We also investigated the maximum sum-rate of all three linear receivers in 
the low SNR regime. When there is feedback, we considered a simple opportunistic scheduling 
algorithm based on exploiting SINR feedback from the receivers to the transmitter. We derived 
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new accurate maximum sum-rate approximations and exact maximum sum-rate scaling laws for 
the MRC, ZF and MMSE receivers. Our results demonstrated that all three linear receivers obeyed 
the same asymptotic scaling law, however the rate of convergence differed based on a simple 
intuitive parameter a, which we quantified. Our results were confirmed through comparison with 
Monte Carlo simulations. 

Appendix 

A. Proof of the SINK C.D.F. with MRC 

The normalized SINR Xmrc can be written as 

|hlfchi,fc| _ Z 

^^^C = |. t . ,2 = V I 1 ^^^^ 

where Z is a chi-squared distribution with 2Nr degrees of freedom with c.d.f. 

Nr-l k 



fe=0 



and y is a chi-squared distribution with 2{Nt — 1) degrees of freedom with p.d.f. [24] 

N 9 

r(A^.-i) (i) 



N 9 

fviy) = ^^W^i ■ (39) 



The c.d.f. of Xmrc can be written as 

^Xmrc (^) = (^MRC < X) 

POO 

= / Friz <xiy + l))fYiy)dy 
Jo 



r(jv, - 1) U\ feo p.« W h 



Finally, we obtain the desired result by first solving the integral in (|40l) using identities in [18] 
followed by some algebraic manipulation. 

B. Proof of the MRC Maximum Sum-Rate with Round Robin Scheduler 
Substituting ([TOl) into ([17]) . we can write the maximum sum-rate of MRC as 

Nr-l k „oo x^k 



RmrM = Yl z]"P'fc / 7 wv^^^ 

k=0 p=0 (l + ^x^ 
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E (!)(-!)'-' --d. (41) 

k=0 p=0 g=0 ^^-^ 

where the first line follows using integration by parts, and the integral in the last line is solved 
using integral identities in [18]. 



C. Proof of MMSE Maximum Sum-Rate at High SNR for Nf > Nr 
We first write the maximum sum-rate in (fT9l ) as 

Nr-l k ^Nt~l\ k 

i?MMSE(p) = E E J(-l)'"'^(/') (42) 

fc=0 p=max(0,fc-Aft+l) 9=0 

where 

e{p) = (^^) V (^1 _ AT, + g, (43) 

and K = Nt — q — l+ p. To obtain high SNR maximum sum-rate expressions for (l42l) . we require 
the limiting behavior of F ^1 — A^^ + q, ^ j as p ^ oo. This is given by [15] 

Hn, r f 1 - A', + ^) ^ (^(iV. - ,) - log 



p— >oo 



1-Aft+g 



A^t-l-? VP/ 

where </)(■) is the digamma function. Substituting (l44l) into (1431) . it is convenient to write ^(p) as 

%) = ^i(p) + ^2(p) + of-l.^ (45) 



where 



and 



''.(P)H'f)"(^^^^(w-,)-log(^|| (46) 



We proceed by first analyzing the high SNR performance of and Q2 in (|46]) and (147]) . For ^1, 
as g < A't — 1 we are only required to find the minimum exponent of in (l46l) . i.e. minimum 
K. This occurs when g = A^,. — 1, so that = Nt — + p. Clearly, for any value of p > 0, 

= o (!^) . (48) 
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We now consider 62{p). For large SNR, we only need to consider the smallest exponent of p in 
(|47l) . i.e. p = 0. This gives 

02{p) = ^r^ T + o(-] . (49) 



Nt-q-1 \p 

Substituting (|48l) and (|49l ) into (|43l) . and then using the resulting expression in (l42l) . we obtain 
the maximum sum-rate as 

'Nt-l\J^ fk\ (-1)^-9 /logp\ ^ 1 /logp 

' 



fc=0 ^ ^ g=0 t >i \ / ^ ' 



(50) 



where the last equality follows after some simple algebraic manipulation. 



D. Proof of Low SNR Results 
From [16], ^ and So are given respectively by 

^^0 min 

# and g,^_21n2[i^-(0)]- 

iVomin i?"(0) i?-(0) 

where R{-) and /?(•) denote the first and second derivative of R{-) respectively w.r.t. p. 
We start with the derivation of ^ and So for the MRC receiver. Substituting (l3]) into (fTSl) 

J»o min 

and taking the derivative w.r.t. p, we find that 

^MRc(O) = log2(e)E [|hj , h,fc|l = Nr hg,{e) . (52) 



Similarly, taking the second derivative of -Rmrc w.r.t. p, we obtain 



^MRCIUJ 



IVi ' "h I 



log2(e)iV,(2iV, + Ar,-l) 



(53) 



where we have used (|38l) and (l39l) . Substituting (|52|) and (l53l) into (ISTI) . we obtain the desired 
result. 

Now consider the MMSE receiver. In this case, it is convenient to write 

^mmse(p) = N,\og, {^-Y. {--w) ^ Hfc (KKt)^h,,j (54) 
which is obtained by substituting ^ into (fT5l) . and using [25]. In this form, we can calculate 
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the required first and second order derivatives as follows 



^mmse(O) = log2(e)E 



KkKk 



and 



To solve (l58l). we write 



E 



h^KKth,, 



Tr (E 



hj , KK^h 



i,k 



K^h,,,ht K 



2iNt - l)Nr 



where we have used [26, Theorem 2.3.5]. Substituting (1571) into (1561) gives 



-^mmse(O) 



log2(e) 



(2(iVi - l)iV, + iV,(iV, + 1)) 



log2(e)iV,(2iV, + iV,-l) 



Finally, substituting (l55l) and (|58l) into (ISTI) . we obtain the desired result. 



(55) 



(56) 



(57) 



(58) 



E. Proof of Theorem [7] 

We first give the following lemma: 
Lemma 2: For sufficiently large K, 

i^-"(p)=log.(l + ^F-(^)) . (59) 
Proof: The proof follows by applying a general result from order statistics [19] to Cx^^{p) 
in (|27]) . and performing some simple algebraic manipulation. ■ 
Lemma [2] implies that for large K, we need only consider the c.d.f. in (flOl) . (|7]) and ([8]) in the 
high X regime. By consider the expressions for the largest exponents of x in (flOl) . ([7]) and dS]), 
followed by some algebraic manipulation, this is given for MRC, ZF and MMSE respectively 



as 



^-x^Nr-Nt I Nt 



1 - 



Nt~l 



(Nr-iy. 



OiX 



Nr-Nt-1-^ 



(60) 
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Fx^A^) = 1 - + oix''^-'''-') (61) 

and 

/ATX Nt-l 

Fx^MsA^) = 1 - e-^;g^._ix^^-^' l^-^j + 0(0;^^-^-^) . (62) 

Note that these c.d.f. expansions are each of the general form F^{x) = 1 — ae'^x^'""^' + 
0(2:^'^"^'"^). In addition, these expansions approach the actual c.d.f. (without the o(-) term) for 
large x. Now, the factor (7^) (l59l ) can be solved by substituting N{K) = K + 1 into the 
iteration (|3TI) . As K is large, only one iteration is needed, since more iterations would produce 
nested In terms which have negligible effect. This gives for large K 

i?--(p) = iV,log2 (1 + 1^ + HK) + {Nr - Nt)0{\n\nK)) 

(63) 



The proof follows by algebraic manipulation. 
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TABLE I 

Absolute error (%) between exact inverse c.d.f. and analytical values using the MMSE receiver for 



DIFFERENT ITERATIONS, WITH K = 4 AND /5 = 3 dB. 



Antenna Configuration 


1 It. 


2 It. 


3 It. 


4 It. 


5 It. 


6 It. 


Nt=Nr = 2 


10.3845 


0.8346 


0.0685 


0.0056 


0.0005 


0.0001 


Nt=Nr=3 


21.6231 


2.7165 


0.3519 


0.0458 


0.0060 


0.0008 


Nt = iV,. = 4 


32.8893 


5.1821 


0.8443 


0.1383 


0.0227 


0.0037 


TVt = iV, = 5 


43.9760 


8.0287 


1.5140 


0.2877 


0.0547 


0.0104 




SNR, p 



Fig. 1. Maximum sum-rate of MRC, ZF and MMSE using scheduling without feedback. Comparison between analytical and 
Monte Carlo simulated maximum sum-rate for Nt = Nr = 4. 



21 




Receive Antennas, N 

r 

Fig. 2. SNR crossover point where the ZF maximum sum-rate is equal to the MRC maximum sum-rate using scheduling 
without feedback. 
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SNR, YdB 



Fig. 4. Percentage of MMSE maximum sum-rate achieved by MRC using scheduling without feedback. 




SNR, YdB 



Fig. 5. Percentage of MMSE maximum sum-rate achieved by MRC for A'',. — 8 using scheduling without feedback. 
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Fig. 6. Maximum sum-rate of MRC and MMSE using sciieduling witliout feedback. Comparison between analytical and Monte 
Carlo simulated maximum sum-rate for different antenna configurations with Nt > Nr using scheduling without feedback. 




Fig. 7. Maximum sum-rate of MRC, ZF and MMSE using the opportunistic scheduler with SINR feedback. Comparison 
between analytical and Monte Carlo simulated maximum sum-rate for different antenna configurations with p — 3 dB. 
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Fig. 8. Maximum sum-rate difference between MMSE and MRC vs. log In if for Nt = Nr — i and p = 3 dB using tlie 
opportunistic scheduler with SINR feedback. 
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